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Abstract 

We extend our approach to the Casimir effect between absorbing dielectric multi- 
layers [M. S. Tomas, Phys. Rev. A 66, 052103 (2002)] to magnetodielectric systems. 
The resulting expression for the force is used to numerically explore the effect of the 
medium dispersion on the attractive/repulsive force in a metal-magnetodielectric 
system described by the Drude-Lorentz permittivities and permeabilities. 
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1 Introduction 

The attractive theoretical possibility of obtaining the repulsive Casimir force 
between (planar) magnetodielectrics [1] can hardly be realized at micron and 
submicron distances because there is no natural medium with permeability 
different from unity in the relevant frequency range [2]. Therefore, as follows 
from the Lifshitz theory [3], the Casimir force should always be attractive at 
distances exploited in modern experiments and applications [4]. However, it 
has (also) been conjectured that for media with nontrivial magnetic properties 
the contributions from the low-frequency side of the spectrum, where perme- 
ability dominates permittivity, may still lead to the repulsive total force [5]. 
In other words, to obtain the Casimir force properly, it is necessary to ac- 
count for the medium dispersion over the entire spectrum and not only in the 
A ~ d range [6]. In view of important potential applications in the develop- 
ment of micro- and nanoelectromechanical systems [7] and nanotechnology in 
general [4] , it seems therefore worthwhile to consider the effect of the medium 
dispersion on the Casimir force between magnetodielectrics in more detail. 
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The aim of this work is to reconsider the statements made in Refs. [1,2] and 
[5] by analysing the attractive/repulsive Casimir force at zero temperature in 
a very simple metal-magnetodieletric system and accounting for the medium 
dispersion. However, in order to rigorously provide an alternative (and more 
familiar) result for the Casimir force to that given in Ref . [1] , we first extend 
the theory of the Casimir effect between absorbing dielectric multilayers [8,9] 
to magnetodielectric systems. 



2 Theory 



Consider a multilayered system described by permittivity e(r, uj) = e'(r, uj) + 
ie"(r,u) and permeability /i(r,uj) = fj,'(r,u) + i/i"(r,uj) defined in a stepwise 
fashion, as depicted in Fig. 1. 
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Fig. 1. System considered schematically. The dashed line represents the plane in a 
lossless layer where the stress tensor is calculated. 



In order to calculate the relevant component of the Maxwell stress tensor Tj z 
[11] 

Tj. zz = — (E Z D Z — En • Dm + B Z H Z — Bii ■ Hi 
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we decompose the macroscopic field operators into positive frequency and 



negative frequency parts according to 

E(r,t)= dwE(r,cu)e- i<JJt + duE^(r,u)t 
Jo Jo 

With the constitutive relations [10] 



iuit 



(2) 



D(r, uj) = e(r, w)E(r, uj) + 47rPjv(r, uj), 

H(r, uj) = — rB(r, uj) - 47rMjv(r, uj), 

fj,{r, uj) 
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they obey the macroscopic Maxwell equations of the standard form. Here 
Pjv(r,tc>) and Mjv(r, a;) are the electric and magnetic noise polarization oper- 
ators, respectively which obey the commutation rules (in the dyadic form): 

[P N (r, uj), P] v (r', uj')} = ^e"(v, uj)I5(v - v')5(uj - 00'), (4a) 

[M N (v, uj), M] v (r', a/)] = A^^ M(r ~ *)K«> ~ <A (4b) 

where I is the unit dyadic. Therefore, any (annihilation) field operator is re- 
lated to Pat(i", uj) and M_ N (r, uj) via the classical Green function G(r, t';uj) 
satisfying 
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V x 1 V x -e(r,u)%-I- 



V>(r,v) ' c 



G(r, r'; uj) = 4vrK(r - r') (5) 



and the outgoing-wave condition at infinity. As a consequence, all field corre- 
lation functions can be expressed through the Green function in accordance 
with the fluctuation-dissipation theorem [12]. In particular, for the electric- 
field correlation function we find 

E(r, cu)E t (r', cu')) = -^ImG(r, r'; uj)5(uj - uj'), (6) 

1 IT C 2 

and the magnetic-field correlation function is obtained from this expression 
using B(r, u) = (—ic/uo)V x E(r, uj). 

Applying the above results to the jth layer and taking into account that £j{oj) 
and fJ,j(u) are real and that P N (r,u) = and ~M N (r,uj) = in this region, 
for the relevant correlation functions in Eq. (1) we find 

h f°° uj 2 
{E(r,t)T>(T,t)) Te(j) = -J ^H-ImG^r.r^), (7a) 



(B(r,t)H(r,t)> ) = - / du—^lmGf (r, r; uj). (7b) 
tuj it Jo fi>j{v) 

Here Gj(r, r'; uj) is the Green function element for r and r' both in the layer 
j, and 

Gf(r,r»=VxG 3 (r,r';w)xV' (8) 
is the corresponding Green function element for the magnetic field. 

The Casimir force per unit area fj between the stacks bounding the jth layer 
equals to Tj jZZ when omitting its infinite-medium part [13]. Formally, this is 
done by replacing the Green function in Eq. (7) with its scattering part 

Gf (r, r'; uj) = G,(r, r'; uj) - G°(r, r'; uj), (9) 
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where G°(r, r'; lo) is the infinite- medium Green function. Writing G^ c (r, r'; 00) 
simultaneously as 



Gf(r,r>)=/.>)Gf(r,r';u;), (10) 

we have 



(ii) 

where 

kj(u) = ^H^-M^ (12) 

is the wave vector in the layer and G^ (r, r'; a;) = G^ xx (r, r'; 00) +G^ yy (r, r'; a;). 

In terms of G* c (r, r'; a;), the above expression for the force formally coincides 
with the corresponding result for a purely dielectric multilayer [8]. As follows 
from Eq. (5), G^ c (r, r'; u) itself is of the same form as the Green function 
Gj c (r, r'; uo) for a purely dielectric system [14], the only difference being that 
the wave vectors in the layers are now given according to Eq. (12) and the 
Fresnel coefficients of the surrounding stacks are modified because of the dif- 
ferent magnetic properties of the system in the present case. Accordingly, with 
these modifications in mind, we can formally adopt all subsequent results ob- 
tained in Ref. [8]. Thus, with Sj{uj) = = 1, according to Eq. (2.15) of 
Ref. [8] we have for the Casimir force 



+Gf„>,r; W )-Gf i >,r;a,}, 



" ■k k J W>\^, D m (ui,k) 



2vr 2 Jo Jo 3 q ^ s D qj (i£,k) 

where f3j(u, k) = \Ju 2 /c 2 — k 2 , 

D qj (oo, k) = l- r]_(u, k)r] + (u, k)e 2i ^ , (14) 

and rj±(u, k) are the reflection coefficients of the right (+) and left (— ) stack 
for TE (q = s)- and TM (q = p)-polarized waves, respectively. The second 
line in Eq. (13) is obtained in the usual way by converting the integral over 
the real a;-axis to that along the imaginary cj-axis, letting u = i£, introducing 
Kj(£, k) = —i/3j(i£, k) = \Ji 2 /c 2 + k 2 , and noticing the reality of the integrand. 
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3 Discussion: Repulsive force 



For numerical calculations, it is convenient to make transition to polar coor- 
dinates in the second line of Eq. (13). Letting £/c = Kj cos0, k = Kj sin0, and 
then x = 2Kjd, we have (hereafter we omit the index j) 



15 f°° , /""V 2 

7o sin i 



^Cj POO p 

f = M d )^ J Q dxx 3 e~ x j 



v r q _(x,(l))rl(x,<f>) 



q=p,s 



where f\d{d) = ir 2 hc/240d A is the Casimir force in the ideal (attraction) case 
[15] and r±(x,<f>) = r±(i||cos0, ^sin0). For a vacuum-medium interface we 
explicitly have 



e — \ ea cos 2 6 + sin 2 6 

r p (x,<j)) = V ==, r s (x,0) =r p [e^/i], (16) 

e + J e/j, cos 2 + sin 2 



where e and /x are functions of uo = i(cx/2d) cos0. Clearly, owing to this 
dependence of e and /i on d, the integral in Eq. (15) is also a (slowly varying) 
function of the separation between the slabs. 

Equations (15) and (16) are equivalent to Eqs. (4) and (5) quoted by Kenneth 
et al.[l]. In the subsequent lines, these authors proceeded by regarding e and 
/i as some constants (appropriate for the relevant frequency range [5]). In 
the following, we reconsider their analysis adopting the single-oscillator model 
for the electric (e) and magnetic (to) polarization [16]. Each slab is therefore 
described by permittivity and permeability of the Drude-Lorentz type 

{e> ^1 = 1 + {^Pud/Cf 

{2ui r p ll d/ c) 2 + x cos <fi(x cos + 2^ v d/ c) ' 

where up u [y = e, to) measure the coupling strengths between the medium 
and the electromagnetic field, whereas ujt v and j u are resonant frequencies 
of the medium and linewidths of the associated resonances, respectively. It 
is immediately seen that the non-dispersive regime is reached at distances 
ujTvd/c ^> 1, where e and /i acquire their static values {e, /i}(0) = \-\-uPp v j u)\ v . 
For convenient measures of the magnitudes of e and /i, we can therefore take 
the quantities up u /u>Tu- To have a feeling about the importance of the second 
term in the denominator of Eq. (17), we note that the main contribution to 
the integral in Eq. (15) comes from the value of x of about 2. 

We have performed numerical calculations for a metal-magnetodielectric sys- 
tem. For the dielectric function em of the metal we set up e = up, u>T e = 0, 
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and 7 e = 7 in Eq. (17) and for the permeability we set /im = 1 (wp m = 0). To 
simulate 7 of good metals (such as Au, Cu, Al) at zero temperature, we have 
let 7 = 10~ 6 7au(300 K) which, in conjunction with ^7a u (300 K) = 35 meV 
[17], corresponds to the residual relaxation parameter of gold at T = K [18]. 
The metal is considered a reference medium and characteristic frequencies of 
the magnetodielectric are scaled with up. Also, the distance between the slabs 
is measured in units kp 1 = c/ujp. For orientation, we note that in the case of 
Au (hu P = 9.0 eV [17]), kp 1 ~ 22 nm and a similar value is found for other 
good metals. Thus, submicron distances correspond to k P d < 50. 




0.5 1 1.5 2 2.5 3 



Fig. 2. Contour plot of ///id as a function of P e = ojp e /u)Te and P m = ujp m / tuxm 
for kpd = 1, ijjTejujp = 0.7, and uJTm/wp = 0.5. Damping parameters are 
7 = 3.9 x 10 -9 u;p and je/cuTe = lm/^Tm = 10~ 2 (these parameters are kept 
fixed in all simulations). Constant-force lines are plotted in steps of A/ //id = 0.02, 
starting from ///id = —0.06 at the top left corner. 

The dependence of the reduction factor ///id on plasma frequencies u)p v (since 
uj T u are fixed) of the magnetodielectric is illustrated in Fig. 2. Note that for this 
d, the figure actually represents the Casimir force in units /o = iv 2 hckp /240 
(~ 5.67 kN/m 2 for Au). One observes that the transition from the attractive 
to repulsive regime occurs roughly when P m > P e and that the constant-force 
lines at larger P m and P e in the transition region are (nearly) straight lines 
characterized by a (nearly) constant ratio P m /P e . Similar behavior, however, 
in terms of e and fi, is found for the border line between the attractive and 
repulsive regions in Fig. la of Ref. [1]. We also note that our simulation agrees 
with the conclusion that for large e and /i the sign and the magnitude of the 
force depend on the surface impedance Z = \Jfj,/e of the magnetodielectric 
[1]. This can be seen if we rewrite Eq. (17) in terms of the relative quantities 
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Pvi Qu = ^Tv/ujp and d = kpd 



{e,/i} = 1 + 



P^QJ? 



(2Q u d) 2 + x cos (j>(x cos + 2Q v d^ u / ujtv) 



At large P u one can disregard unity here, so that Z ~ P m /P e . Accordingly, 
(nearly) straight lines seen in the (P m , P e ) contour plot correspond to constant- 
Z lines. Since for this separation between the slabs e and n are strongly disper- 
sive (note that kx e d = 0.7 and kT m d = 0.5), one must conclude that accounting 
for the medium dispersion does not much affect the conclusions drawn in Ref. 
[1] regarding the (e, /^-dependence of the force. 




Fig. 3. ///id as a function of t^e/Vp and ujTm/^p for kpd = 1, ujp e /ujT e = 0.5, and 

UPm/uJTm = 3. 

Clearly, of crucial importance for the appearance of the repulsive force is the 
relative position of the electric and magnetic resonance of the magnetodi- 
electric. Figure 3 illustrates the variation of ///id with transverse frequencies 
(Q u ) of the magnetodielectric while keeping the magnitudes of e and \i (deter- 
mined by P v ) fixed. At large Q v the force tends to a constant value because 
{e, fi} — ► 1+P„ in Eq. (18). Of course, in the opposite limit, the force vanishes 
as {e,fi} ~ 1 for small Q v . It is seen that the appearance and magnitude of 
the repulsive force is governed by the position of the magnetic resonance for 
almost all Q e . Note that in the region where Q e > Q m , which is normally 
always the case, the force is attractive at this distance. 

Since Q v and d appear in the form of the product Q u d = u Tu d/c in Eq. (18), 
the mismatch between the electric and magnetic resonances can be compen- 
sated by properly adjusting d and thus one can in principle still observe the 
repulsive force although at larger distances. This possibility is illustrated in 
Fig. 4. One observes that, indeed, even for very small values of Q m (<^C Q e ) at 
larger distances there is a repulsive force between the media. To estimate the 
distance at which the repulsive force could be observed in a (possibly) realistic 
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Fig. 4. ///id as a function of the distance kpd between the slabs and the posi- 
tion uJTm/up of the magnetic resonance for ujp e /u>Te = 0.5, u)Te/ojp = 0.7, and 

m — 

system, we note that the highest frequencies known so far at which perme- 
abilities differ significantly from unity are (for some antiferromagnetics) in the 
THz range [20]. Realistic values of the parameters Q u are therefore Q e ~ 0.1, 
corresponding to ojTe in the range of optical frequencies, and Q m ~ 10 , cor- 
responding to ujTm m the THz range. The distance dependence of the (total) 
force for these parameters across the repulsive region is presented in Fig. 5. 

2 - 
1 - 
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Fig. 5. / in units 10 _14 /o as a function of the distance kpd between the slabs for 
a "realistic" system: w^e/wp = 0.1, uJT m /ojp = 10 -4 and other parameters are the 
same as in Fig. 4. 

As seen, the force is attractive at small distances, which supports claims made 
in Ref [2] . With increasing d, however, the effect of permeability is enforced and 
the force crosses over to the repulsive regime already at kpd ~ 640 (which is 
~ 14/xm for Au). Thus, the effect of magnetic polarization appears at distances 
smaller by two orders of magnitude than it would be expected on the basis 
of a simple estimate d ~ Ay m = 2n x 10 4 kp 1 . Evidently, this effect is due to 
the cumulative contribution of the waves from the low-frequency side of the 
spectrum and is therefore a direct consequence of the medium dispersion at 
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long wavelengths. To estimate the magnitude of the repulsive force at such 
distances, we calculate its largest value occurring at kpd ~ 850, where we 
find / ~ —2.88 x 10~ 14 /o- Accordingly, the maximal repulsive force between, 
e.g. Au and this hypothetical [uo Te = 0.9 eV, oo Pe = 0.45 eV, 7 e = 9 meV, 
^Tm = 0.9 meV, up m = 2.7 meV, and 7 m = 9 /xeV] magnetodielectric is 
/ ~ —0.16 nN/m 2 and is observed at the distance d ~ 18.7 jum. 

Of course, at higher temperatures the above estimate becomes invalid as at 
such large separations between the slabs the thermal corrections must be ac- 
counted for. The force at nonzero temperature is formally obtained from Eq. 
(13) by letting [9] £ -> £ m = 2nmk B T/h and f™d£ -> (2ixk B T/h) £~ =0 (1 - 
5 m o/2), where £ m are the Matsubara frequencies. One may immediately con- 
clude that the force in the present system is always attractive at room tem- 
perature. Indeed, at T = 300 K and for Au, the mth contribution to the force 
is governed by the quantity 

= 2nm^-k P d = 1.8 x 10- 2 m k P d. 
c njjjp 

Thus, at distances kpd > 10 2 only the m = term contributes as all other 
terms are strongly exponentially damped. Since for the Drude metals r s M (0, k) = 
[22], this term depends only on the dielectric properties of the magnetodi- 
electric and the force is therefore attractive. We have verified this by di- 
rectly calculating the force at T = 300 K and, for example, at kpd ~ 850, 
where the zero-temperature force is maximally repulsive, we have obtained 
/ ~ 5.09 x 10- 13 /o ~ 2.89 nN/m 2 . 

The validity of the above estimates rests upon the adequacy of the Drude 
(local) model for the metal adopted. This model is often used in the consid- 
erations and precise calculations of the Casimir force between metals and is 
generally believed to adequately describe the properties of good metals at low 
frequencies [17,21]. Strictly speaking, however, the electromagnetic response 
of a metal is nonlocal, at least in the frequency range of the anomalous skin 
effect, and should therefore be described by the nonlocal reflection coefficients 
r q M [11]. The effect of the metal nonlocality below uj < uj p on the Casimir 
force between two metal (Au) plates at zero temperature has recently been 
considered in detail in Ref. [23] and it has been found that anomalous skin 
effect produces only a minor correction (below 0.5% for d above 50 nm) to 
the force calculated using the local description of the metal (which also agrees 
with the assertion given in Ref. [21]). Moreover, as can be concluded from 
Fig. 9 of Ref. [23], this correction diminishes and saturates with increasing 
distance between the plates. Accordingly, accounting for the anomalous skin 
effect of the metal cannot significantly change the above estimates on the re- 
pulsive force at T = K. As follows from the small frequency behaviour of 
the surface impedances [23], since the property r s M (0, k) = holds nonlocally 
as well, the force at T = 300 K is, as before, always attractive. However, since 
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also r^f(0,k) 7^ 1, contrary to the generally accepted r^(0,k) = 1, further 
consideration of the effect of the metal nonlocality on the Casimir force at 
finite temperatures is needed. 

We end this short discussion by noting that very recently Henkel and Joulain 
reported a similar analysis [24] , however, paying more attention to the depen- 
dence of the force on geometrical rather than on material parameters as we 
have done. In this respect, these two works are complementary. As concerns 
our findings on the distance dependence of the force, these are (generally) in 
agreement with those of Henkel and Joulain concerning thick media. 



4 Summary 

Using the properties of the macroscopic field operators appropriate for mag- 
netodielectric dissipative systems, we have extended our previous approach to 
the Casimir effect between two dielectric multilayered stacks to magnetodi- 
electric systems. As expected, the resulting expression for the force between 
the stacks is of the same form as for a purely dielectric system. However, 
the Fresnel coefficients are modified owing to the different magnetic proper- 
ties of the system in the present case. Using this result, we have numerically 
analysed the effect of the medium dispersion on the Casimir force in a metal- 
magnetodielectric system described by the Drude-Lorentz permittivities and 
permeabilities. According to this analysis, taking into account the medium 
dispersion does not much affect the conclusions that one may draw from a 
simple nondispersive model of Ref. [1] regarding the material dependence of 
the force. Our simulations also demonstrate that the force at zero temperature 
in a realistic system is attractive at small distances (covering the submicron 
and micron range for good metals), which supports statements in Ref. [2], and, 
depending on the magnetic (dispersion) properties of the magnetodielectric, 
becomes repulsive at larger distances. However, owing to contributions of the 
low-frequency waves to the force, the crossover occurs at a much smaller dis- 
tance than expected from the simple d ~ Xxm estimate. At room temperature, 
the force on good (Drude) metals such as Au, Cu, Al remains attractive over 
the entire range of distances. 
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